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1 Introduction 



The Calogero-Moser-Sutherland (CMS) and the Ruijsenaars-Schneider-van Diejen (RSvD) 
interacting many-particle models play a distinguished role among the integrable Hamiltonian 
systems, having numerous relationships to important fields of mathematics and physics. They 
have profound applications in the theory of solitons (see e.g. pQ, [2], [3], [1]), and recently they 
appeared in the context of random matrix theory as well (see e.g. [5], [0]). Quite surprisingly, 
these intriguing relationships are well-understood only for the models associated with the An 
root system. It appears that the main technical obstacle for developing analogous theories in 
association with the non-A^-type root systems is the lack of knowledge of explicit action-angle 
variables for the non-A„-type CMS and RSvD models. In this paper we wish to narrow this 
gap by constructing action-angle systems of canonical coordinates for both the hyperbolic BCn 
Sutherland and the rational BCn RSvD models. 

In order to define the i?C„-type hyperbolic Sutherland and the rational RSvD many-particle 
systems, we first introduce the subset 

c = {x = (xi, . . . , x„) G I xi > . . . > Xn > 0} C M", (1.1) 

which can be seen as an appropriate model for the open Weyl chamber of type BCn- Let us 
recall that the phase space of the Sutherland model is the cotangent bundle of c, for which we 
have the natural identification 

= cxW = {{q,p)\qe c,peW}. (1.2) 

Recall also that the hyperbolic BCn Sutherland dynamics is generated by the interacting many- 
particle Hamiltonian 



^'^tp?-^ E (- 

c=l l<a<b<n ^ 



,sinh {qa - qb) sinh (g^ + qh) 

(1.3) 

n 2 " 2 

+ V - + V - 

^sinh2(g,) ^sinh2(2g,)' 

where the so-called coupling parameters g, gi and g2 are arbitrary real numbers satisfying the 
inequalities g^ > and gf + g2 > 0. In other words, we are interested only in the hyperbolic 
BCn Sutherland mode0 with purely repulsive interaction. Since the strength of the interaction 
is governed by the numbers g"^, g\ and g"^-, they are usually called the coupling constants. 

Though the non-A^-type CMS models have received a lot of attention in the last couple of 
decades (see e.g. the fundamental papers [7j, [8] and the book |9]), the symplectic reduction 
understanding of the BCn Sutherland model with three independent coupling constants is a 
quite recent development [lOj. Besides providing a Lax representation of the dynamics, the 
symplectic reduction approach has the advantage that it naturally leads to a fairly simple 



^Notice that with the speciahzation g2 — ^ we recover the hyperbohc -B„ Sutherland particle system, mean- 
while with gi — we obtain the C„-type model. 
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solution algorithm of purely algebraic nature. By pushing forward the reduction picture, in 
this paper we are able to furnish action-angle variables for the standard BCn Sutherland model 
with repulsive interaction. 

The non-A^-type deformations of the classical Ruijsenaars-Schneider many-particle systems 
have been introduced by van Diejen [TT]. Just as for the Sutherland model, the phase space of 
the rational BCn RSvD model is the cotangent bundle T*c, which is naturally identified with 
the manifold 

= cxW = {{X,9)\Xe c,9 eW}. (1.4) 
Recall that the rational RSvD dynamics is characterized by the Hamiltonian 



/ 2\f/ 2\f"/ 



(Ac-A,)V V (Ac + A,)^ 

('i^c) (1.5) 



VilV' ■ A^y v 



where fi, v and k are arbitrary real parameters satisfying yU 7^ 7^ z/. Although the Liouville 
integrability of the non-A^-type RSvD models has been verified pJL], the Lax representation 
of their dynamics is still missing. However, by generalizing our results on the C„-type RSvD 
model [12], in this paper we provide a Lax matrix and an elementary solution algorithm for the 
rational RSvD model (11. 5p with z/k > 0. Moreover, the proposed reduction approach permits 
us to construct action-angle variables as well. 

The organization of the paper can be outlined as follows. Section |2] is devoted to a brief 
account on the necessary group theoretic and symplectic geometric background underlying the 
derivation of the Sutherland and the RSvD models from a unified symplectic reduction frame- 
work. In Section [3] we review the symplectic reduction understanding of the BCn Sutherland 
model. Although this is a standard material (see [lO]), our new contribution on the spectral 
properties of the Lax matrix of the Sutherland model, formulated in Lemma [H seems to be 
crucial in advancing the reduction approach to cover the RSvD model, too. Starting with Sec- 
tion H] we present our new results on the rational BCn RSvD model. By fitting the BCn RSvD 
model into a convenient symplectic reduction picture, we are able to provide a Lax matrix and 
an elementary solution algorithm as well. The main technical result is Theorem [5l in which 
we confirm that the parametrization of the Lax matrix of the RSvD model does provide a 
Darboux system on the reduced phase space. In Section |5] we elaborate on the consequences 
of the proposed reduction approach. In particular, a natural construction of canonical action- 
angle variables for both the Sutherland and the RSvD model comes for free. Furthermore, the 
action-angle duality between the repulsive BCn Sutherland model and the BCn RSvD system 
with z/K > becomes also transparent. 

This paper is a continuation of our recent work [12] on the hyperbolic Sutherland and the 
rational RSvD models associated with the C„ root system. It is a very fortunate situation 
that many results for the models associated with the BCn root system can be derived almost 
effortlessly by generalizing the analogous results of the C^-type particle systems. Since in [12] 
we have carried out a very detailed analysis on the particle systems of type C„, in this paper 
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we can be brief on many aspects of the i?C„-type models. Though our presentation tries to 
be self-contained, in this paper we rather focus on differences between the C.„-type and the 
i?C„-type models, and we provide proofs only for those facts that have no natural analogs in 
the Cn case. Therefore the reader may find it useful to have a copy of [12] on hand while 
reading the paper. 



2 Preliminaries 

In this section we gather the necessary group theoretic and symplectic geometric material 
underlying the unified symplectic reduction derivation of the hyperbolic BCn Sutherland and 
the rational BCn RSvD models. Throughout the paper our group theoretic conventions try to 
be consistent with the book [13J, whereas the symplectic geometric conventions come mainly 
from IMj. To facilitate the comparison with our work on the C„-type models, the majority of 
the notations are directly borrowed from paper [T2] . 

Take an arbitrary positive integer nGN={l,2, ...} and let = 2n. With the aid of the 
N X N unitary matrix 

1 

In 



c 



E U{N) (2.1 



we define the non-compact real reductive matrix Lie group 

G = U{n, n) = {ye GL{N, C) | y*Cy = C}. (2.2 
The corresponding real matrix Lie algebra has the form 

g = u(n, n) = {Y e Ql{N, C) | Y*C + CY = 0}, (2.3 

on which the map 

( , ) : g X ^ M, (Fi, Fa) ^ (Fi, F2) = ^(^1^2) (2.4 

provides a symmetric Ad-invariant non-degenerate bilinear form. 

Let us remember that the fixed-point set of the Cartan involution Q{y) = (y*)^^ {y E G 
naturally selects a maximal compact subgroup 

K = {yEG\ e{y) = y} = {yEG\y is unitary} = U{n) x U{n) (2.5 

of the Lie group G. Also, the Lie algebra involution 6{Y) = —Y* {Y E q) corresponding to O 
naturally induces the Cartan decomposition 

g = «©p (2.6 

with the Lie subalgebra and the complementary subspace 

« = ker(^ - Id) = {F G g I r* = -Y} and p = ker(^ + Id) = {F G g | F* = F}, (2.7 

respectively. Due to the Cartan decomposition (12. 6p . each F G g can be decomposed uniquely 
as 

F = F+ + F_ {Y+Et,Y_E^). (2. 
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Next, notice that the set of diagonal matrices 

a={Q = diag(gi, . . . , -gi, . . . , -g„) G p | (gi, . . . , g„) G R"} (2.9) 

forms a maximal Abelian subspace in p (12. 7p . Let a"*" denote the subspace of the off-diagonal 
elements of p; then we have the orthogonal decomposition p = a © a"*". Let us also consider the 
centralizer of a inside K, which is the Abelian subgroup 

M = ZK{a) = {diag(e'^\ . . . , e^^", e'^\ . . . , e'^") G | (xi, • • • , Xn) G M"}. (2.10) 

Obviously its Lie algebra has the form 

m = {diag(ixi,. . . , ixn, ixi, • • • , iXn) e ! | (xi, • • ■,Xn) e R"}. (2.11) 

If m-*- denotes the subspace of the off-diagonal elements of t, then we can write t = m © m^. 
Finally, recalling the Cartan decomposition (12. 6p . we end up with the refined decomposition 

= m©m-^ © o© o-^. (2.12) 

Having equipped with the above group theoretic objects, in the rest of the section we review 
some basic notions from symplectic geometry. Recall that the cotangent bundle T*G of the Lie 
group G can be trivialized, say, by left translations. Upon identifying the dual space q* with 
the Lie algebra g via the bilinear form (12. 4p . it is clear that the product manifold 

V = GxQ = {{y,Y)\yeG,Ye5} (2.13) 

provides a convenient model for T*G. Furthermore, the tangent spaces of V can be naturally 
identified as 

T^y,Y)V = T(^y,Y){GxQ)^TyG®TYQ = TyG®g {{y,Y) eV). (2.14) 

Let us observe that on the model space V = T*G the canonical one-form G Q^(V) reads 

{6y © 6Y) = {y-'6y, Y) {{y, Y) eV,6y®SYe TyG © g), (2.15) 

whereas for the canonical symplectic form we use the convention u = — d?? G Q^{V). 

Now note that the smooth left action of the product Lie group K x K on the group manifold 
G defined by the formula 

{kL,kR).y = kLyk^' {yeG,{kL,kn)eKxK) (2.16) 

naturally lifts onto T*G. Working with the model space V (12.131) of the cotangent bundle, the 
lift of the above K x K-action (I2.16P takes the form 

[kL, kn) . {y, Y) = {kLyk-^\ kRYk-^^) ((y, Y) G (fc^, kn) eKxK). (2.17) 

This action is clearly symplectic, admitting the K x i^'-equivariant momentum map 

J: P ^ (€©«)*=«© (y, F) ^ J(y, Y) = {yYy-^)+ © (-r+). (2.18) 
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Without any further notice, in the rest of the paper we shall frequently use the natural dual 
space identification © t)* = £ © 6 induced by the bilinear form fl2.4p . 

To proceed further, with each column vector V G subject to the conditions V*V = N 
and CV + V = Owe associate the Lie algebra element 

^{V) = ifi{VV* - In) + i(/i - i^)C e t, (2.19) 

where /i, u are arbitrary real parameters satisfying fi ^ ^ u. Also, let E G denote the 
distinguished column vector with components 

Ea = -En+a = l (aGN„ = {l,...,n}), (2.20) 

and consider the Lie algebra element 

Jo = (-^(E)) ©i/tC G «©e, (2.21) 

where k is an arbitrary real parameter. In order to derive the hyperbolic BCn Sutherland and 
the rational BCn RSvD models from symplectic reduction, we wish to reduce the symplectic 
manifold {V, u) at the very special value Jq 02.211) of the momentum map J fl2.18p . We mention 
in passing that the parametrizations of the Lie algebra elements ^(V^) fl2.19p and Jq f l2.2ip turn 
out to be very natural in the sense that, after performing the reduction, the parameter triple 
(/i, u, k) can be identified with the coupling parameters of the rational BCn RSvD model (11. 5p . 

Our experience with the CMS and the RSvD models convinces us that, in general, the 
application of the shifting trick leads to a shorter and neater derivation of these particle systems 
from symplectic reduction. As the initial step of the shifting trick (see e.g. [15]), we have to 
identify the adjoint orbit passing through — Jo (I2.2ip . Since C commutes with each element of 
K, for the adjoint orbit in question we have the natural identification O © {—inC} = O, where 

O = 0{^{E)) = {^{V) et\V eC^, V*V = N,CV + V = 0}. (2.22) 

Following the prescription of the shifting trick, we also introduce the extended phase space 

p-t = VxO = {iy,Y,p) \yeC,Y eQ,peO}, (2.23) 

and endow it with the product symplectic structure 

w"'^* = u + u^, (2.24) 

where, of course, w'^ is the standard Kirillov-Kostant-Souriau symplectic form carried by the 
orbit O f l2.22p . The natural extension of the K x i^-action fl2.17p onto V^^^ is given by the 
diagonal action 

{kL,kR) . iy,Y,p) = ikLyk],\knYk],\kLpkl'), (2.25) 
and the corresponding K x i^-equivariant momentum map takes the form 

jct ^ pcxt ^^^^^ Y, p) ^ r'^'iy, y, p) = {{yYy-')+ + p) © {-Y^ - ikC). (2.26) 
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As a matter of fact, it is clear that J'^^* takes its values in the subalgebra 

s{t®i) = {XL®XRet®t\ tr(Xi) + tr(X^) = 0} < « © (2.27) 

Now, the shifting trick guarantees that 

V//j,{K xK)= V'^y/oiK X K), (2.28) 

i.e. for our purposes it is an equally valid approach to perform the Marsden-Weinstein reduction 
of the symplectic manifold ("P*^^*, cj'^^*) at the zero value of the momentum map J"^^^ (12.261) . 

3 The hyperbolic BCn Sutherland model 

In this section we review the Lax matrix and the symplectic reduction understanding of the 
hyperbolic BCn Sutherland model with three independent coupling constants. Our discussion 
on the reduction aspects of the model is mainly based on the ideas presented in [TD], adapted 
to the conventions of [12]. 

3.1 The Lax matrix of the Sutherland model 

The main goal of this subsection is to introduce the Lax matrix of the Sutherland model, and 
to analyze some of its spectral properties that prove to be pertinent in the symplectic geometric 
understanding of the BCn RSvD model. As a preparatory step, with each point {q,p) G we 
associate the diagonal matrices 

Q = diag(gi, . . . , g„, -qi, • • • , -Qn) e a and P = diag(pi, . . . -pi, -pn) G a. (3.1) 

Let adg denote the restriction of the linear operator adg = [Q, ■] G 0t(0) onto the off-diagonal 
part of the Lie algebra q (12.31) . Notice that the regularity condition g G c ensures the invertibility 
of the linear operator adg. Therefore, making use of the standard functional calculus, the matrix 

Lp{q,p) =P- sinh(a~dQ)-^e(^) + coth(adQ)(i/€C) G p (3.2) 



is well-defined. Let us note that the above introduced N x N matrix Lp = Lp{q,p) is Hermitian 
with block-matrix structure 

~ A B' 
-B -A ' 



(3.3) 



where A and B are n x n matrices satisfying A* = A and B* = —B. More concretely, for their 
matrix entries we have 

. — iu . „ iu ■ \ 

Aa,b = ^-TT T, ^c,c = Pc, Ba,b = ^-TJ ] 7, ^c,c = ■ , s + lKCOth(2gJ, 

smh(ga - qb> smh(ga + qb) smh(2gc) 

(3.4) 

where a,b,c & N„ and a ^ b. Now, with the parameter triple [fi, u, k) we associate the map 

L:V' ^5, (g, p) ^ L{q, p) = Lp (g, p) - inC. (3.5) 
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As one can see in [TO], the above map L provides a Lax matrix for the hyperbohc -BC„ Suther- 
land model. The exact relationship between the parameters (/i, z/, n) and the Sutherland cou- 
pling parameters (9,91,92) appearing in f ll.3p will be clarified later (see f l3.65p and fl3.66p ). 

Having defined the Lax matrix of the BCn Sutherland model, we now turn our attention 
to its spectral properties. Remembering that the only difference between L and Lp is the 
ant i- Hermit ian constant term inC, it is clear that the spectral properties of the non-Hermitian 
Lax matrix L fl3.5p can be understood by analyzing the spectrum of the Hermitian matrix Lp 
(13. 2p . Since Lp belongs to the complementary subspace p (12. 7p . we know from general principles 
that it can be conjugated into the maximal Abelian subspace a (12. 9p by some element of the 
maximal compact subgroup K (12. 5p . However, this diagonalization procedure becomes much 
more explicit by exploiting the singular value decomposition of the sum of the matrices A and 
B introduced in the block- matrix decomposition (13. 3p . More precisely, we can write 

A + B = usv*, (3.6) 

where u and v are n x n unitary matrices, meanwhile s = diag(si, . . . , s„) is a diagonal matrix 
filled in with the singular values si > . . . > s„ > of the matrix A + B. Now, upon defining 
the N X N block-matrices 



"'2 



V + U V — u 

V — U V + u 



and S 



s 
-s 



(3.7) 



one can easily verify that k ^ K, S ^ a and Lp = kSk ^. Having diagonalized the matrix Lp 
(13. 2p . from the definition (13. 5p we see at once that 

L^ = k{S^ - KHN)k-\ (3.8) 

therefore the spectrum of the Hermitian matrix L^ can be identified as 

a{L^) = a{S^ - kHn) = cr{s^ - kH^). (3.9) 

On the other hand, remembering the singular value decomposition (13. 6p . we can also write 

{A + B){A- B) = {A + B){A + B)* = usv*vsu* = us\-\ (3.10) 

from where we obtain the spectral identification 

a{s^) =a{A^ - B^ -[A,B]). (3.11) 

Now, the comparison of the equations (13.91) and (13. lip immediately leads to the formula 

a{L^) = a{A^ - B^ - [A, B] - kH^). (3.12) 

Since is an Hermitian matrix, the spectral mapping theorem guarantees that each eigenvalue 
of L is either a real number or a purely imaginary number. However, under certain technical 
assumptions, the relationship (I3.12p permits us to provide a more accurate description for the 
spectrum of L. 
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Lemma 1. Suppose that z/ 7^ 2/i and > 0; then for each point {q,p) G we have L{q,pY > 
0, i.e. the matrix L{q.,pY is positive definite. In particular, the eigenvalues of the Lax matrix 
L{q,p) are non-zero real numbers. 

Proof. Take an arbitrary point {q,p) G and keep it fixed. First, notice that if k = 0, tlien 
tfie Lax matrix L = L{q,p) is of type C„. However, we liave a fairly complete knowledge on the 
spectrum of the Lax matrix of the C„ Sutherland model. Namely, since i> 7^ 2/i, from Lemma 
1 in we see that the Hermitian matrix L is invertible, whence > is immediate. 

In the following we assume that k 7^ 0. As an important auxiliary object in our proof, let 
us consider the -BC„-type Lax matrix L = L{q,p) associated with the parameters (/i, u — n, 0). 
Recalling (13.51) and (13. 3p . we see that L is an Hermitian matrix with block- matrix decomposition 



A B 
-B -A 



G p. (3.13) 



Moreover, remembering (13.41) . for the matrix entries of A and B we have 

sinh(g<j - g;,) ' "'^ sinh(g„ + g;,) ' sinh(2gc)' 

where a, 6, c G N„ and a ^ h. Since L is Hermitian, it is clear that U' > 0. Thus the direct 
application of (I3.12p on the Lax matrix L yields immediately that 

- B^ -[A,B]>0. (3.15) 

In order to find the connection between the Lax matrix L (13.51) and the auxiliary Lax matrix L 
(13.131) . we introduce the diagonal matrix q = diag(gi, . . . ,g„). Comparing the equations (13. 4 p 
and (13.141) . it is obvious that 

A = A and B = B + in coth{q) . (3.16) 

Thus it is immediate that 

A'^-B'^ - [A,B] -k^In 

= A'^-B^- [A, B] + ifi:[coth(qr). A] - i/t(coth(qr)5 + B coth(q)) + smh{q)-^. (3.17) 

Upon introducing the column vector V G with components Vc = 1/ sinh(gc) (c G N„), the 
right hand side of the above equation can be simplified considerably. Indeed, by applying the 
standard hyperbolic identity 

coth(x) + coth(y) ^ 1 

sinh(a; + y) sinh(x) sinh(j/) ' 

one can easily verify the relations 

iK{coth{q)A - A coth(q)) = -fii^VV* + fiK sinh(qr)-^ (3.19) 
iK{coth{q)B + B coth(q)) = -finVV* + (^/t + k^) sinh(q)-l (3.20) 
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Plugging these formulae into fl3.17p . we end up with the concise expression 

A^-B^- [A, B] - k^In = A^-B^- [A, B] + z//€sinh(q)~2. (3.21) 

However, due to f lS.lSp and the assumption uk, > 0, the matrix on the right hand side of the 
above equation is manifestly positive definite. Therefore, according to the spectral identification 
flXT2]) . we conclude that > 0. □ 

3.2 The phase space of the Sutherland model 

In this subsection we perform the Marsden-Weinstein reduction of the extended symplectic 
manifold (P*^^*, w*^^*) at the zero value of the momentum map J^^^ fl2.26p . As a first step of the 
reduction, we have to solve the constraint 

J^^*(y,r,p) = (3.22) 

for {y, Y, p) G V"^^^. In other words, we have to understand the differential geometric properties 
of the closed level set 

£o = (^'^"*)"'({0}) = {iy,Y,p) e V"^' I r''\y,Y,p) = 0} C P'^^*. (3.23) 

In order to derive the phase space of the Sutherland model from the proposed reduction picture, 
we are looking for a special parametrization of £9 induced by the KAK decomposition of the 
group elements y & G. (For background information on the KAK decomposition see e.g. the 
book p^.) As it can be seen from the lemma below, besides the KAK decomposition, the 
most important ingredient of the parametrization is the Lax operator (13. 5p . 

Lemma 2. Suppose that u + k, ^ 0; then for each point {y, Y, p) G -Cq of the level set there are 
some g G c, p G M" and rii.rjji G K , such that 

y = VLe'^VR\ Y = 7]RL{q,p)r]^\ p = r]L^{E)r]l\ (3.24) 

Proof. Take an arbitrary point {y,Y,p) G -Cq- The KAK decomposition tells us precisely that 
the Lie group element y E G can be decomposed as 

y = kLe'^k],\ (3.25) 

where kiykji E K and 

Q = diag(gi, . . . , g„, -gi, . . . , -g„) G a (3.26) 

with some gi > . . . > > 0. Also, by fl2.22p . we can write p = C,{V) with some column vector 
V eC^ satisfying V*V = N and CV + V = 0. 

Plugging the above parametrizations into the constraint J'^^^{y,Y,p) = 0, the explicit form 
of the momentum map J^^^ fl2.26p immediately leads to the relationship y+ = —inG, together 
with 

= iyYy-')+ + p = kL (sinh(adQ)(A;^V_M - cosh(adQ)(ifi:C) + ^kl'V)) kl\ (3.27) 
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Upon introducing the shorthand notations 

Y_ = k^^Y_kR G p and V = k^^V G C^, (3.28) 
from the equation fl3.27p it readily follows that 

sinh(adQ)F_ = -^{V) + cosh(adQ)(i/«C) = -^{V) + iKcosh{2Q)C. (3.29) 
Spelling out the components of the above matrix equation, for all k, I G Na^ we have 

sinh(gfc - qi){Y_)k,i = -ifx{VkVi - 5k,i) - iifi - i^)Ck,i + incos\i{2qk)Ck,u (3.30) 

where it is understood that qn+c = —qc for all c G N„. 

Now take an arbitrary c G N^. With the specialization k = I = c the above equation fl3.30p 
takes the form 

o = -iMIK|'-i), (3.31) 

whence it is obvious that = ~^+c = e'-^" with some real parameter Xc £ IR- Notice also that 
with k = c and I = n + c the equation fl3.30p translates into 

sinh(2gc) (^-)c,n+c = iz^ + iKcosh(2gc)- (3.32) 

Therefore, under the assumption u + k, ^ 0, from the above relationship it is clear that q^ 7^ 0. 

Next, let a, 6 G N„ be arbitrary numbers satisfying a ^ b. With the specialization k = a 
and / = b the equation fl3.30p has the form 

sinh(g, - qb){Y_)a,b = -ifiVaVb ^ 0, (3.33) 

therefore qa 7^ qb- Putting the above considerations together, we see that > . . . > g„ > 0, 
whence the regularity property g G c is immediate. 
Now let us introduce the diagonal matrix 

m = diag(e^^S . . . , e^^", e'^S . . . , e'^") G M. (3.34) 

Due to the construction of m, it is clear that m~^V = E (12. 20 p . Therefore, by applying the 
linear operator Ad^-i on the equation (13. 29 p . we obtain 

sinh(adQ)Ad™-i(yi) = -^{E) + cosh(adQ)(iKC). (3.35) 

By inspecting the diagonal and the off-diagonal parts of Adm,-i(y_) separately, and utilizing 
the regularity of q, we can write 

Ad^-i(Y'_) = P- sinh(adQ)-i^(E) + coth(adQ)(i/tC) (3.36) 

with some diagonal matrix P = diag(pi, . . . —pi, . . . , —pn) G ci. Thus, remembering the 
definition (13. 2p . we can simply write = mLp{q,p)m~^; therefore the relationship 

F_ = kRmLp{q,p)m~^k^ (3.37) 

is immediate. Since V = k^V = kimE, we also have 

p = ^{V) = kLmi{E)m~^kl\ (3.38) 

Therefore, with the group elements til = kim G K and rjR = kum G K, the lemma follows. □ 
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To proceed further, let us notice that the Abehan group 

f/(l), = {(e^^^l^, e'n^) eK xK\xeM} = f/(l) (3.39) 

is a closed normal subgroup of the product Lie group K x K, whence the coset space {K x 
K)/U{1)^: inherits a natural (real) Lie group structure from K x K. Let us also consider the 
smooth product manifold 

= V^ X {K X K)/U{1),. (3.40) 

Having equipped with the above objects, now we can introduce a natural parametrization of 
the level set £o (13.231) motivated by Lemma [2l Indeed, by imitating the proof of Lemma 2 in 
[12], one can easily verify that the map 

T^: ^p-t, {q,p,{r^L,VR)U{l).) ^ (vlc'^Vr ,VRL{q,p)VR ,VLaE)r]j:') (3.41) 

is a well-defined injective immersion with image T'^(A^'^) = £o- Also, just as in the proof of 
Lemma 3 in [12], the parametrization T"^ makes it easy to verify directly that the zero element 
of the Lie algebra s{t(Bi) (I2.27P is a regular value of the momentum map J^^* (12. 26 p . Therefore 
the level set 2o is an embedded submanifold of V^^^ in a natural manner. More precisely, there 
is a unique smooth manifold structure on £,q such that the pair {£,q, lq) with the tautological 
injection 

Lo-.Ho^V^''' (3.42) 

is an embedded submanifold of At this point let us notice that, due to the relationship 

T'^(A^'^) = lq{£,o), the map T'^ factors through (£o, '•o); therefore there is a well-defined map 

T^:7W^^£o, (3.43) 

such that T"^ = ° ^o- Let us observe that, since lq is an embedding, the map Tq is automat- 
ically smooth (see e.g. Theorem 1.32 in [E]). Now, since Tq is a smooth bijective immersion 
from Ai^ onto and since it acts between manifolds of the same dimension, it is immediate 
that T^ is a diffeomorphism. The above ideas can be summarized by saying that the diagram 




(3.44) 



is commutative. In other words, the pair {Ai^ , T'^) provides an equivalent model for the smooth 
embedded submanifold {2o,lo). 

Utilizing the model {Ai^ , T^) of the level set in the following we complete the symplectic 
reduction of (P''^*, o;*^^*) at the zero value of the momentum map J*'^*. For this purpose let us 
note that on the model space Ai^ (I3.40p the residual K x i^- action takes the form 

(kL^kR) .{q,p,{r]L,VR)Uil)*) = {q,p, {kiVL, kRr]R)U{l)^). (3.45) 
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Therefore it is obvious that the orbit space Ai^/{K x K) can be naturally identified with the 
base manifold of the trivial principal {K x ir)/f/(l)^,-bundle 

T,s. j^s ^ p5 (^^^^ {r]L,VR)Uil).) ^ (3.46) 

An immediate consequence of the above observation is that the reduced symplectic manifold 
can be identified as 

V'^y/oiK xK)= M^/{K xK)= . (3.47) 

As is known from the theory of symplectic reductions, the reduced symplectic form G Vt^{V^) 
is uniquely determined by the condition 

{-K^fuj^ = {T^Yu"'''. (3.48) 

However, since the derivative of ( KAT\ can be worked out explicitly, the computation of the 
pull-backs in f l3.48p is almost trivial. Either doing the calculations by hand, or remembering 
the results presented in [TU], the following theorem is immediate. 

Theorem 3. The reduced symplectic form can be written as = 2 Yl^=i ^^Ic ^ ^Pc- That is to 
say, up to some trivial rescaling, the globally defined coordinate functions qc, Pc (c G N^) form 
a Darboux system on the reduced manifold . 



3.3 Solution algorithm for the Sutherland model 

The goal of this subsection is to present a solution algorithm for a class of Hamiltonian systems 
in association with the family of the Ad-invariant smooth functions defined on the Lie algebra 
Q. For, take an arbitrary Ad-invariant smooth function F : g — ?■ M, i.e. we require 

F{yYy-^) = F{Y) (VF G 0, G G). (3.49) 

Now, let 

prg:P""* = G'xgxO^0 (3.50) 

denote the canonical projection onto g; then it is obvious that the composite function pr*F = 
Fopr^ is a smooth K x i^'-invariant function on p*^^*; therefore it survives the reduction. More 
precisely, the corresponding reduced Hamiltonian (pr^F)*^ G C°°{V^) has the form 

(pr*F)^ = FoL, (3.51) 

as can be readily seen from the defining relationship 

(n'riwlFf = {T'TwlF. (3.52) 

It is a standard fact in reduction theory that the Hamiltonian flows of the 'unreduced' Hamilto- 
nian system (P*^^*, u"^^^, pr*F) staying on the level space £o project onto the flows of the reduced 
system {V^,u^,F o L). However, flnding the 'unreduced' flows is a relatively simple exercise; 
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therefore this projection method gives rise to a natural and efficient solution algorithm for the 
reduced Hamiltonian system. 

To make the above observation precise, we need the integral curves of the Hamiltonian vector 
field Xpr*F G X{V^^^) generated by pr*F. Recalling w*^^* fl2.24p . from the defining relationship 

Xpr*^ja;^^* = d(pr;F) (3.53) 

we find easily that at each point {y, Y, p) G P^^* the Hamiltonian vector field has the form 

{X^r;F)(y,Y,p) = (z/VF(r)), © Oy © 0, G T(,,y,,)P^^\ (3.54) 

where the gradient VF(Y) G g is defined by the condition 

(VF(F), 6Y) = (dF)y(5F) (V5r G TyQ = fl). (3.55) 

From fl3.54p it is clear that the induced Hamiltonian flows are complete, having the form 

R3t^ (yoe*^^(^»\ 1^0, Po) G V"^' (3.56) 

with some {yQ,YQ,po) G V'^^^. Therefore the reduced Hamiltonian flows are also complete. 
Now take an arbitrary flow 

R3t^ {q{t),p{t)) eV^ (3.57) 

induced by the reduced Hamiltonian F o L, and let Lq = L(g(0),p(0)). Since the unreduced 
flow passing through the point {e'^^^\ Lq,^{E)) E £,q at t = projects onto the reduced flow 
(I3.57p . from the definition of the parametrization T*^ fl3.4ip it is clear that for each t G M there 
is some {r]L(t),r]R(t)) E K x K such that 

^gQ(o)gtvF(Lo)^ Lo, m) = Mt)e'^^SRitr\ Mt)mt),p{t))Mt)-\ r^L{mE)vL{t)-'). 

(3.58) 

By comparing the 0-components of the above equation we see that 

Lo = r]R{t)L{q{t),p{t))r]R{t)-\ (3.59) 

from where we conclude that during the time evolution of the reduced dynamics the Lax matrix 
L (13. 5 p undergoes an isospectral deformation. What is even more important, the G-component 
of equation f l3.58p leads to the relationship 

gQ(o)giVF(Lo)eiVF(Lo)*eQ(o) ^ r]L{t)e"^'-'^r]L{t)-\ (3.60) 
which entails the spectral identification 

Somewhat more informally one can say that the matrix flow Q{t), and so the trajectory q{t), 
can be recovered simply by diagonalizing the matrix flow 
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Besides providing a nice solution algorithm, the above observation can also be seen as the start- 
ing point of the scattering theoretic analysis of the reduced Hamiltonian system {V'^ , u^, FoL). 
Indeed, by analyzing the temporal asymptotics of the matrix flow fl3.62p . one can understand 
the temporal asymptotics of the trajectory q{t) as well. 

In order to establish the connection with the Sutherland many-particle systems, notice that 
the Hamiltonian of the BCn Sutherland model (11.31) with three independent coupling constants 
can be realized as the reduced Hamiltonian induced by the Ad-invariant quadratic function 

F,{Y) = \{Y,Y) = \tT{Y') (Yeg). (3.63) 

Indeed, a simple computation reveals that 

(pr*F2)^ = i/^ (3.64) 



with coupling constants 



Note that if vk, > 0, then the interaction is purely repulsive. Conversely, if one starts with 
an arbitrary triple of the non-negative Sutherland coupling constants {g'^,gf,g2) satisfying the 
inequalities > and gf + g2 > 0, then the corresponding repulsive Sutherland model can be 
recovered e.g. by choosing the parameter triple (/i, u, k) with components 

I I 1^721 + ^/gl + Wl '^^gl 

V2 |(72| + \/gl + %f 

Notice that the above defined parameters satisfy the inequalities /i < 0, > and k > 0, 
therefore the assumptions made in Lemmas [1] and [2] are automatically met. As for the solution 
algorithm of the Sutherland model, note that VF2(Lo) = whence from equation (I3.62p 

we see at once that the trajectories can be determined by diagonalizing the matrix flow 

t^^Q^^)e¥^^^e¥^l. (3.67) 



We close this section with some remarks on the range of the attainable Sutherland coupling 
constants coming from the proposed symplectic reduction picture. From equation ( I3.66P it 
is clear that any non-negative triple of the coupling constants isP' ■, g\i g^) with g^ > ^ and 
g\ ^ g\ > can be realized by an appropriate choice of the parameters (/x, z/, h) satisfying 
VK > 0. Due to the repulsive nature of the interaction, in these cases the Sutherland model has 
only scattering states. It conforms with Lemma [1], according which the eigenvalues of the Lax 
operator L are real. Thus the scattering properties of the model can be understood explicitly 
by analyzing the temporal asymptotics of the truly exponential-type matrix flow (13.671) . 

However, looking at (13. 65 p . it is obvious that from the proposed reduction picture we can 
derive the Sutherland model even with gl < 0, provided that uk < 0. Since Lemma [1] does not 
apply in these cases, in certain non-empty region of the phase space some of the eigenvalues 
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of the Lax operator may be purely imaginary. Therefore, as can be conjectured from f l3.67p . 
some of the Sutherland particles may exhibit oscillatory behavior with bounded trajectories. To 
fully understand the details of this phenomenon, one should sharpen Lemma [T] to characterize 
the spectral properties of the Lax matrix for uk, < 0. Nevertheless, since our primary interest is 
to study the duality properties of the standard hyperbolic Sutherland BCn model with purely 
repulsive interaction, we leave this interesting exercise for a future study. 



4 The rational BCn RSvD model 

In this section we work out the rational BCn RSvD model with three independent coupling 
constants from the proposed symplectic reduction framework. Let us keep in mind that Section 
[3] and the present section provide the necessary technical background to establish the action- 
angle duality between the standard repulsive Sutherland model (11. 3p and the rational RSvD 
(II. 5p model with uk > 0. Therefore, according to our discussion in the concluding remarks of 
Sectional we may assume at the outset that fi < 0, u > and /t > 0. 



4.1 The phase space of the RSvD model 

First of all, we need some objects introduced in the study [12] of the C^-type rational RSvD 
model associated with the pair of non-zero parameters (/i, u). For each a G N„ let us consider 
the complex-valued rational function 

. 3 A ^ ( V . . . . A„) . ..(A) ^ - (l . ^) n (l . (l + ^) e C, ,4,1) 

Recall also that the Lax matrix of the rational C„ RSvD model is defined by the matrix-valued 
function 

A:V^^exp{p), {X,e) ^ AiX,e), (4.2) 
where the matrix entries lying in the diagonal n x n blocks are given by the formulae 

Aa,biX,e) = e'^^'''\za{\)z,i\)\'^ , (4.3) 

ZljJ + Aa — Ab 

A JX P] - r-^a-db ^aW^bW 2i/x 

\Za{A)Zb{A)\'2 ^1/^ - Aa + Ab 

meanwhile the matrix entries belonging to the off-diagonal n x n blocks have the form 

Aa,n+t{X, 0) = = e'-~'^Zb{\)M\)zb{\V\^ ^. ^'1" , , + '^^^^Kb. (4.5) 

2l/i + Aa + Afe l/i + Aa 

for any a, 6 G N„. Besides the Lax matrix A^ from the theory of the C„ RSvD model we need 
the column vector J-'(A, 6) G with components 



Fa{\e) = e''^\za{\)\-- and 7n+a{\0) = e-"" Za{X)\za{\)\-\ (4.6) 
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where a G N„. Also, in the forthcoming computations we shall frequently encounter the column 
vector 

V{\,e) =A{\,9)--^T{\,9) eC^. (4.7) 

Notice that V(A,^) is well-defined, since the positive definite Lax matrix A{X,9) G exp(p) has 
a unique square root belonging to exp(p). Furthermore, as it can be seen from Proposition 8 
in [ig, we have V*V = and CV + V = 0. 

In order to handle the BCn RSvD model associated with the parameter triple (//, z/, k), we 
need some new objects, too. In particular, let us introduce the smooth functions a and (3 
defined on the positive half-line by the formulae 



aix) = ■;= and nix) = ik—= — , (4.j 



where x G (0, oo). Now, with each positive n-tuple A = (Ai, . . . , A^) G (0, oo)" we associate the 
n X n diagonal matrix A = diag(Ai, . . . , A^) and consider the Hermitian N x N matrix 



h{X) 



a{\) /3(A)' 
-/3(A) a(A) 



(4.9) 



Making use of the functional equation a{xY + /3(x)^ = 1, one can show that h{X)Ch{X) = C, 
i.e. the matrix h{X) is an Hermitian element of the Lie group G (12. 2p . Finally, let us introduce 
the shorthand notation 

A = diag(Ai, . . . , A„, -Ai, . . . , -A„) G a. (4.10) 

Having equipped with the above objects, we are now in a position to provide an appropriate 
parametrization of the level set £,q fl3.23p based on the diagonalization of the Lie algebra part 
of £o. 

Lemma 4. Suppose that v ^ 2^, z/ + k 7^ and vk > 0; then for each point {y, Y, p) G £0 there 
are some A G c, G and rji^rjn G K , such that 

y = VLA{X,e)h{Xr'r^^\ Y = r^nh{X)Ah{X)-'vR\ p = VLmKO))vl' . (4.11) 



Proof. Take an arbitrary point {y,Y,p) G £o- Remembering the momentum map J°^^ (12.2 
it is clear that 1^+ = —inC. On the other hand, since any element of the subspace p (12. 7p can 
be conjugated into a (12.91) by some element of K (12. 5p . we can write 

r_ = kuDk],' (4.12) 

with some kn & K and 

D = diag((ii, . . . , (i„, -di, . . . , -dn) G a (4.13) 
satisfying di > . . . > dn > 0. Therefore 

Y = kniD - iKC)k^\ (4.14) 
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from where we obtain 

= kniD^ - KH^)k],\ (4.15) 

However, by combining Lemmas [1] and [2], it is obvious that Y"^ is positive definite, i.e. Y"^ > 0. 
Giving a glance at the above equation fl4.15p . it is thus clear that = a/c?^ — is a well- 
defined positive real number for each a G N„, satisfying the inequalities Ai > . . . > A„ > 0. 
Now, recalling h{X) f l4.9p and utilizing the functional equations 



a(x)'' - /3(xf = W 1 + ^ and 2a(x)/3(x) = — , (4.16) 
\ x^ X 

it is not hard to see that 

h{\)Ah{\)-^ = D- inC. (4.17) 
Plugging this observation into fl4.14p . we get 

Y = kRh{X)Ah{X)-^k^\ (4.18) 

To proceed further, notice that the momentum map constraint yields the relationship 

= {yYy )+ + p = + p (4.19) 

as well. This is clearly equivalent to the equation 

y*yY-Y*y*y + 2y*py = 0. (4.20) 

However, due to (12.221) we can write p = ^(V^), where V G is an appropriate column vector 
satisfying V*V = N and CV + V = 0. Thus, recalling ( 12.190 . the above constraint (14.201) can 
be cast into the form 

y*yY - Y*y*y + 2ipiy*VV*y - 2ifiy*y + 2i(^ - u)C = 0. (4.21) 

Plugging the parametrization fl4.18p into the above equation, we get 

2ifxhk*j^y*ykRh + Ahk*j^y*ykRh - hky*ykRhA = 2ifx{hk*^y*V){hky*Vy + 2i(/i - u)C. (4.22) 

At this point let us note that the last equation can be naturally identified with equation (31) 
in [16]. Since this equation was the starting point of our scattering theoretic analysis of the 
hyperbolic C„ Sutherland model, we have full control over the structure of its ingredients. In 
particular, due to Lemma 1 in [16|, we know that A must be a regular element of a, whence 
Ai > . . . > A„ > 0, i.e. A G c. Furthermore, by Lemma 2 in [IB], we can write 



h{X)ky*ykRh{X) = mA{X, 9)m* (4.23) 

with some 6 G MJ^ and m G M. Noticing that the group element h{X) G G commutes with each 
element of M fl2.10p . we obtain 

{ykRmh)*ykRmh = m*hk*Ry*ykRhm = A. (4.24) 
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Third, our analysis yields the relationship 

h{X)ky*V = mJ^{X, 6) (4.25) 

as well. Now let us notice that equation fl4.24p allows us to give a characterization of the global 
Cartan decomposition (polar decomposition) of ykRmh{\). Indeed, we have 

ykRmh{\)=TiLA{\d)-- (4.26) 

with some rji ^ K; thus the parametrization 

y = r]LA{\,e)hi\y\kRmy' (4.27) 

is immediate. Plugging this into fl4.25p and remembering the definition f l4.7p . we obtain at once 
that V = ?7iV(A,6'); therefore the relationship 

P = aV)=VLm>^,0))vl' (4.28) 

also follows. Finally, upon setting rjR = kum G K, from the equations (14.270 . ( I4.18P and fl4.28p 
we see that the proof is complete. □ 

Motivated by the above lemma, we can introduce a natural parametrization of the level set 
£o (I3.23p . To make this idea precise, first let us define the smooth product manifold 

= V^ X {K X K)/U{1),. (4.29) 

Now, a trivial generalization of the proof of Lemma 10 in [12] immediately convinces us that 
the smooth map 

defined by the assignment 

(A, 6, ivL, VR)Um ^ {vlA{X, e)h{Xr'r^],\r^nh{X)Ah{X)-'r^],\ VLm\ 0))vl') (4.31) 

is a well-defined injective immersion with image T^{M^) = Ho- Repeating the same arguments 
that we applied in the Sutherland picture, it is clear that factors through the embedded 
submanifold {Hq, Lq), and the resulting smooth map T^: Ai^ — )• £0 is a diffeomorphism. That is 
to say, the pair {Ai^, T^) provides an equivalent model for the smooth embedded submanifold 
('2o, i-o). 

In order to complete the reduction of V^^^ at the zero value of the momentum map J'^^*, let 
us observe that the residual K x i^-action on the model space Ai^ fl4.29p takes the form 

(fci, kn) . (A, 6, {r]L, m)U{l).) = (A, 6, {kL^L, Mi?)f/(1)*). (4-32) 

Therefore the orbit space Ai^/ [K x K) gets naturally identified with the base manifold of the 
trivial principal {K x K)/f/(l)^.-bundle 

T.R-.M''^ P^, (A, e, {r]L, Vr)U{1).) ^ (A, 9). (4.33) 



19 



It is thus evident that for the reduced symplectic manifold we have the alternative identification 

P""7/o(i^ xK)= M^/{K xK)= V^. (4.34) 

Our remaining task is to make explicit the reduced symplectic form G Q'^(V^) naturally 
induced on the reduced phase space V^. Just as in the Sutherland case, it is a tempting idea to 
compute the reduced symplectic structure via a formula analogous to f l3.48p . However, due to 
the presence of the square root of A in f l4.3ip . this approach seems to be hopeless. Instead, 
it is more expedient to invoke the alternative machinery presented in Subsection 4.3 of paper 
|12] . Namely, by analyzing the Poisson brackets of the auxiliary K x 7^- invariant functions 
defined in equations (4.65) and (4.66) of |T2], an almost verbatim computation as in the C„ 
case convinces us that the standard coordinates of are canonical. 

Theorem 5. Utilizing the global coordinate functions Ac and 6c (c G N„) defined on the reduced 
manifold V^, the reduced symplectic structure takes the form = 2 Yl^=i A dX^. 

4.2 Solution algorithm for the RSvD model 

In this subsection we work out an efficient solution algorithm for the rational -BC„ RSvD model 
built on the projection method naturally offered by the symplectic reduction framework. For, 
take an arbitrary K x i^-invariant real- valued smooth function / : G — > M defined on the Lie 
group G, i.e. we assume that 

fiikL, M . y) = fik^yk],') = f{y) (Wy G G, W^k^, k^) e K x K). (4.35) 

Let 

pr^:P'=-t = Gx0xC^G' (4.36) 

denote the canonical projection onto G; then the composite function pr^/ = / o pr^ is clearly 
K X i^-invariant on the unreduced phase space "P^^*. Thus, based on the defining formula 

(vr^)*(pr^/)^=(T^)>^/, (4.37) 

it is easy to see that the corresponding reduced Hamiltonian (pr^/)^ G C°°{V^) has the form 

(pr^/)«(A, 9) = f{Ai\, e)h{\)-') ((A, 9) G V^). (4.38) 

As we have discussed it in the Sutherland picture, the essence of the projection method is that 
the fiows of the unreduced Hamiltonian system staying on the level space £o project onto the 
fiows of the reduced Hamiltonian system. However, from the defining formula 

Xpr^/^a;<=^* = d(pr^/), (4.39) 

we see immediately that at each {y, Y, p) G V^^^ the Hamiltonian vector field ^pr^/ generated 
by pr^/ has the form 

(Xp,^^)(,,y,,) = 0, © (-V/(y))y © Op G r(,,y,,)P^^\ (4.40) 
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where the gradient V/(y) G g is defined by the requirement 

{Vf{y),y-'6y) = {df)y{6y) {\/6y G T,G). (4.41) 

Therefore the Hamihonian fiows of pr^/ are complete, having the very simple form 

R3t^ {yo,Yo - tV/(yo),Po) G (4.42) 

where (|/o,lo)Po) ^ 'P'^^^ is an arbitrary point. Hence the reduced fiows are complete as well. 
Now take an arbitrary fiow 

R3t^ {X{t),e{t)) eV^ (4.43) 

of the reduced Hamiltonian system (V^,u^, (pr^/)^). For simphcity let us now introduce the 
shorthand notations 

A = ^(A(0),^(0)), ho = h{\{0)), Ao = A(0), Vo = V(A(0),^^(0)). (4.44) 
It is obvious that the unreduced fiow 

RBt^ {Alh^\hoAoh,'-tVf{Alh,'),aVo)) e £0 (4.45) 

projects onto (14. 43 p . Recalling T-^ fl4.3ip . we see that for each t G M we can find some pair of 
group elements (?7l(^), '7i?,(^)) E K x K such that 

Alh^' = VL{t)A{X{t), e{t))h{X{t))-'r^n{t)-\ (4.46) 

hoAoK' - tVfiAlh,') = 7^nit)hiXit))Ait)hiXit))-'vRitr\ (4.47) 

e(Vo) = VLimivim, omvLitr'- (4.48) 

From (14.461) we conclude that 

h^'Aoh,' = r]R{t)h{Xit))-'A{X{t), 9{t))h{Xit))-'r]R{t)-\ (4.49) 

which entails the spectral identification 

a(/io Mo/i^') = (r{h{X{t))-^A{X{t),9{t))h{X{t))-^). (4.50) 

Thus, during the time evolution of the reduced Hamiltonian system, the positive definite Her- 
mitian matrix 

A^^iX, 6) = h{X)-^A{X, e)h{X)-^ G G (4.51) 
undergoes an isospectral deformation. Meanwhile, from fl4.47p it follows that 

a{hoAoh,' - tVfiAlh,')) = a{A{t)) = {X,{t), . . . , A„(t), -X,{t), . . . , -A„(t)}, (4.52) 
whence the trajectory 1 1— )■ X{t) can be recovered simply by diagonalizing the linear matrix fiow 

t ^ hoAoho' - tV/(^|^o (4.53) 
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It is worth mentioning that, due to its hnearity in t, the temporal asymptotics of the above 
matrix flow (I4.53P can be analyzed by elementary perturbation theoretic techniques. Note that 
this observation could serve as the starting point of the scattering theoretic analysis of the 
reduced Hamiltonian system (V^jU^, (pr^/)^). 

Our considerations so far apply to any reduced system associated with a K x i^-invariant 
smooth function / f l4.35p . Note, however, that under the assumption uk > the rational BCn 
RSvD model with three independent coupling constants can be nicely fitted into this picture. 
Indeed, upon introducing the K x i^-invariant function 

fi{y) = lti{yy*) (y E G), (4.54) 

one can verify that the corresponding reduced Hamiltonian coincides with the RSvD Hamilto- 
nian (11. 5p associated with the coupling parameters (/i, z/, k), i.e. 

(pr^/i)« = H^. (4.55) 

Let us observe that the assumption uk, > automatically guarantees the lower bound > n 
on the RSvD Hamiltonian. We mention in passing that the verification of (I4.55P is a quite te- 
dious, but elementary calculation. Nevertheless, it can be done easily by utilizing the functional 
identities collected in the appendix of fn\ . 

Turning to the solution algorithm based on (I4.53p . notice that for the gradient (I4.4ip of the 
function fi we have 

V/i(y) = I {y*y - {y*yV) g p {ye G). (4.56) 

Therefore, from (I4.5ip and (14.531) we see at once that the trajectories of the rational BGn RSvD 
model can be determined by diagonalizing the matrix flow 

t ^ hoAoK' - \t [Af" - {Af^r^) , (4.57) 

where A^^ = ^^^(A(0), 6^(0)). We find it remarkable that the properties of dynamics generated 
by the highly non-trivial Hamiltonian (II. 5p can be captured by analyzing the linear matrix 
flow fOTj) . 

To sum up, we see that under the assumption z/k > the rational BGn RSvD model (II. 5p 
can be derived from an appropriate symplectic reduction framework. Parallel to our discussion 
in Section [3] on the Sutherland model (II. 3p with gf < 0, we expect that the RSvD model with 
UK, < can also be understood from symplectic reduction by generalizing Lemma [TJ We wish 
to come back to this issue in a later publication. 

5 Discussion 

In the previous two sections we derived both the standard hyperbolic BGn Sutherland and the 
rational BGn RSvD models from a unified symplectic reduction framework. The derivation of 
the Sutherland model relies on the KAK decomposition of the Lie group part of the level set £o 



22 



(13.231) . meanwhile the symplectic geometric understanding of the RSvD model builds upon the 
diagonalization of the Lie algebra part of £,q. Thereby, by performing the Marsden-Weinstein 
reduction of the symplectic manifold [V^^^ , u'^^^) at the zero value of the momentum map J*^^* 
(I2.26p . we end up with two equivalent realizations, and V^, of the same symplectic quotient 
V°^^//o{K X K). Thus it is obvious that there is a natural symplectomorphism S: — )■ 
making the diagram 



'-0 



'-0 



M^^^S.o^^M'' (5.1) 



R 

p5 = ^jyR 

commutative. In the rest of this section we examine some of the immediate consequences of 
the dual reduction picture (15. ip . First, let us define the functions 

Ae = 5*Ae, = s*e,, = {s-'yq,, = {s^yp,, (5.2) 

where c G N„. Since iS is a symplectomorphism, from Theorems |3] and [5] it follows that 



S*cu^ = 2^d9cAd\c and = (5"^)*a;^ = 2 ^ dgc A dp^- (5.3) 



c=l c=l 



In other words, the globally defined functions Ac and 6c provide a new Darboux system on the 
Sutherland phase space , meanwhile the global coordinates Qc and pc are canonical on the 
RSvD phase space V^. Now we show that these new families of canonical coordinates give rise 
to natural action-angle variables for the Sutherland and the RSvD models, respectively. 

Starting with the Sutherland side of the dual reduction picture, take an arbitrary point 
(g,p) G and let (A,^) = S{q,p) G V^. Now, recalling the parametrizations T"^ (I3.4ip and 
(I4.3ip . from the commutativity of the diagram (15. ip it is clear that 

(e«,L(g,p),e(i^;)) = (r7L^(A,^^)^MA)-^^\r/«MA)AMA)-V>^Le(V(A,^))r^Z') (5-4) 

with some group elements ?7l,77_r G K. By inspecting the g-component of the above equation 
we obtain the spectral identification 

a{L{q,p)) = (t(A) = {±Xc{q,p) \ c G NJ. (5.5) 

Next, take an arbitrary real-valued Ad-invariant smooth function F: g — t- M (I3.49P defined 
on g, and consider the naturally associated reduced Hamiltonian system (V^ , co^, (pr*F)"^). Due 
to (I3.5ip and (15. 4p it is clear that 

(pr;F)^ = F o L = F o A = F(diag(Ai, . . . , A„, -Ai, . . . , -A^), (5.6) 
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i.e. the reduced Hamiltonian depends only on the coordinates Ac (c G N„). It follows that the 
global coordinates Ac and dc (c G N„) provide canonical action-angle variables for the reduced 
system (P'^, w^, (prgF)"^). Note that the action and the angle coordinates of the Sutherland 
picture are exactly the pull-backs of the canonical positions and the canonical momenta of the 
Ruijsenaars picture. 

To conclude the study of the Sutherland side of the dual reduction picture (15. ip . let us recall 
that the hyperbolic BCn Sutherland model can be realized as the reduced Hamiltonian system 
generated by the quadratic Ad- invariant function F2 f l3.63p . Therefore the above construction 
of action-angle coordinates applies to the repulsive Sutherland model equally well. Also, from 
(15. 5p we see that the positive eigenvalues of L (13. 5p provide n functionally independent first 
integrals in involution. Furthermore, the matrix L naturally enters the solution algorithm of 
the Sutherland model (see equation (I3.67P ). whence the non-Hermitian matrix L is indeed a 
Lax matrix for the Sutherland many-particle system. 

In the following we turn our attention to the Ruijsenaars side of the dual reduction picture 
(15. ip . For, take an arbitrary point (A,^) G and let (g,p) = iS~^(A,6') G . Recalling the 
mappings (I4.3ip and T"^ (13.411) . it is clear that 

{A{\,9y^h{\)-\h{\)^h{X)-\i{y{\9))) = {,^Le%-^\^RL{q,p)rR\vLmrL') (5.7) 

with some group elements rj^^riR G K. Now, remembering the definition of the positive definite 
matrix A^'^ (I4.5ip . notice that the G-component of the above equation immediately leads to 
the relationship 

^^C(A, 6) = hixyUix, e)hix)-' = vrc'^^vr, (s.s) 

from where we get the spectral identification 

aiA^'^iX, 6)) = {e^^-/^^^'^) | c G N,}. (5.9) 

We see that ^^""(A,^) has a simple spectrum, and the positive eigenvalues of the Hermitian 
matrix ln(^^^)/2 are exactly the coordinate functions Qc (c G N^). 

To proceed further, take an arbitrary K x i^-invariant smooth function /: G — i- M (I4.35p . 
and consider the naturally generated reduced Hamiltonian system {V^, u^, (pr^/)"^). Recalling 
(I4.38P and (15. 7p . for the reduced Hamiltonian we have 

(pr^/)^(A, 9) = f{A{X, 9)h{X)-') = /(diag(e«^(^'^), . . . , e^"(^'^), e'^^^^'^), . . . , e'^-^^'^))), 

(5.10) 

i.e. the reduced Hamiltonian depends only on the coordinates (jc (c G M„). Thus, it is immediate 
that the global canonical coordinates Qc and Pc (c G N„) form an action-angle system for the 
mechanical system {V^, u^, (pr^/)'^). Let us observe that the action and the angle coordinates 
of the Ruijsenaars picture are coming from the pull-backs of the canonical positions and the 
canonical momenta of the Sutherland picture. 

Now remember that the reduced Hamiltonian system {V^,u^, [pr^fi)^) generated by the 
K X i^- invariant function fi (I4.54p coincides with the rational BCn RSvD model (II. 5p with three 
independent coupling constants. Therefore the canonical coordinates qc and pc (c G N„) provide 
action-angle variables for the RSvD model as well. Notice also that the positive definite matrix 
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A^'^ fl4.5ip plays a distinguished role in the theory of the rational RSvD model. Indeed, the 
positive eigenvalues of the Hermitian matrix ln(^^*-')/2 give rise to n functionally independent 
first integrals in involution. In particular, the matrix A^'^ undergoes an isospectral deformation 
during the time evolution of the RSvD dynamics. Also, remember that A^'^ naturally appears 
in the solution algorithm of the model, as can be seen in equation f l4.57p . Therefore A^'^ fl4.5ip 
meets all the criteria to call it the Lax matrix of the rational BCn RSvD model. 

To sum up, we constructed action-angle systems of canonical coordinates for both the re- 
pulsive hyperbolic BCn Sutherland and the rational BCn RSvD models with three independent 
coupling constants. The relationships between the coupling parameters of the corresponding 
particle systems are displayed in equations fl3.65p and f l3.66p . As we have seen, the action and 
the angle coordinates of the Sutherland model can be naturally identified with the canonical 
positions and the canonical momenta of the RSvD model, and vice versa. That is to say, mak- 
ing use of the dual reduction picture flS.ip . we established the action-angle duality between the 
standard -BC„-type Sutherland and RSvD models. This interesting phenomenon was originally 
discovered by Ruijsenaars in the context of the An-tjpe particle systems [18]. Using advanced 
techniques from symplectic geometry, in the last years the An-tjpe dualities have been reinter- 
preted in the reduction framework, too (see the papers [11], [20]). It appears to be an attractive 
research problem for the future to generalize these techniques to the non-A„-type setup. 

We conclude the paper with some remarks on the possible applications of our results. Besides 
the natural appearance of the Sutherland and the RSvD many-particle systems in the soliton 
scattering description of certain integrable field theories (see e.g. [1], [2] [3], [1]), we expect that 
our results find applications in the theory of random matrices as well. Indeed, by exploiting 
the existing dualities between the A^-type particle systems, the authors of the recent papers [S] 
and [S] have introduced new classes of random matrix ensembles with novel spectral statistical 
properties. Built on the Lax matrices L (13. 5p and A^'^ (I4.5ip . the proposed dual reduction 
picture (15. ip seems to be indispensable in initiating the study of the integrable random matrix 
ensembles associated with non-A„-type root systems. 
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